Charge carrier localisation in disordered graphene nanoribbons 
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Abstract. We study the electronic properties of actual-size graphene nanoribbons subjected to 
substitutional disorder particularly with regard to the experimentally observed metal-insulator 
transition. Calculating the local, mean and typical density of states, as well as the time-evolution 
of the particle density we comment on a possible disorder-induced localisation of charge carriers 
at and close to the Dirac point within a percolation transition scenario. 

^ ' Introduction 

cd; 

Transport in graphene nanoribbons (GNRs) is strongly affected by disorder effects which can 
f^j . be traced back, for example, to dislocations or charged impurities in the substrate, to adatoms 
adsorbed at the graphene surface, to edge defects, or to ripples associated with the soft struc- 
ture of graphene. When modelling disorder, many theoretical studies resort to the generic 
Anderson model [1], which exhibits a disorder-induced localisation transition in three dimen- 
sions (3D) that is absent in lower dimensions, however. One-parameter scaling theory predicts 
that all states are localised for the infinite Anderson-disordered 2D system [2]. The recently 
Q ! observed metal-insulator transition in hydrogenated graphene [3], disordered GNRs [4], and Si- 
MOSFET inversion layers [5] is beyond reach of the Anderson model, but might be explainable 
by percolation-based approaches [6, 4]. 

The classical (geometric) problem of percolation consists in finding a connected path of 
CN ■ accessible sites that spans the whole lattice. On the honeycomb lattice, for site-percolation, this 
will be the case for a concentration of accessible sites p > p c « 0.697 [7]. For real materials the 
strict distinction between accessible and blocked sites seems to be too simplistic. For example, 
^ \ upon hydrogenation, the 7r-bonds of some carbon atoms within a graphene sheet will be blocked 
c3 ■ just partly. Also the electron-hole puddles resulting from charged impurities in the substrate 
lead to a finite difference between the on-site potentials only. Then, tunneling effects between 
the puddles may become possible, allowing for transport despite the absence of a percolating 
cluster. In the quantum case, a spanning cluster does not guarantee transport since scattering at 
its irregular boundaries causes interference effects which may lead to localisation of the charge 
carriers. 

To analyse the localisation properties of low- dimensional systems within numerical ap- 
proaches, both sophisticated algorithms and highly efficient implementations are mandatory 
since the relevant length scales are exceptionally large. In this regard the finite extension of 
mesoscopic graphene flakes or GNRs deserves further attention since it may mask localisation 
effects if the localisation length exceeds the system size. 

In this work we investigate the electronic properties of GNR eigenstates in the vicinity 
of the Dirac point using a local distribution approach based on exact diagonalisation (ED) 
techniques [8]. This method has proven its reliability for the study of Anderson localisation and 
quantum percolation on various lattices in different dimensions. [9, 10] 



Model and method 



To model disordered GNRs we consider the tight-binding Hamiltonian 

N 

H = -* E( c k- + H - c -) + E e A c J (!) 

(ij) i=i 

on a honeycomb lattice with N sites, including electron transfer t between nearest neighbours 
(ij) only. In Eq. (1), c] (q) creates (annihilates) an electron at lattice site i. Drawing the on-site 
potentials from the bimodal distribution 

p[ej] =pS(ej - e A ) + (1 -p)8(e j - e B ) , (2) 

sites are occupied by an atom of type A [B] with probability p [(1 — p)] (binary alloy analogy). 
Without loss of generality we choose the on-site energy of the majority sub-band 6 a = 0, use 
A = 6b — hereafter, and let i fix the unit of energy. 

Localisation properties of disordered systems can be discussed in terms of the local density 
of states (LDOS), 

N 

Pi (E) = J2mn)\ 2 S(E-E n ), (3) 

71=1 

where \i) = c}|0), and \n) is a single-electron eigenstate of H with energy E n . The LDOS can 
be determined very efficiently by the Kernel Polynomial Method which is an expansion of the 
rescaled Hamiltonian into a finite series of Chebyshev polynomials [11]. Thereby an energy level 
broadening appears that can be controlled by the expansion order. 

Within the local distribution approach one has to analyse the behaviour of the normalised 
LDOS distribution f[pi/p me \ upon finite-size scaling for many realisations of disorder (here 
Pme = (Pi) is the mean DOS). While extended states are characterised by a system-size inde- 
pendent f[pi/p me ] : the distribution for localised states strongly depends on TV; its maximum 
shifts towards zero and finally f[pi/p me ] becomes singular as TV —> oc. Since in many cases 
f[Pi/Pme] closely resembles a log-normal distribution [ ], a simplified discussion may use the 
typical DOS, p ty = e^ n<yP ^\ which is directly related to the maximum of the log-normal distri- 
bution. 

Alternatively we may access the localisation properties of a system by determining the 
recurrence probability Pr(£ —> oc), which in the thermodynamic limit is finite for localised states 
and scales to zero as 1/N for extended states. Again a finite series expansion into Chebyshev 
polynomials is promising, in this case applied to the time evolution operator [9]. Then, for 
example, we may track how an initially localised wave packet evolves in time by calculating the 
time dependent local particle density, 

I N 2 

l^r^^l^l^e-^H^O))^^) . (4) 

m=l 



Numerical results 



As compared to mesoscopic graphene flakes, the DOS of ordered GNRs exhibits a series of spikes 
(van Hove singularities) indicating quasi- ID behaviour. The number and position of these spikes 
depend on the ribbon geometry and on the (finite) number of unit cells in the transverse cross 
section. Furthermore, GNRs form very interesting edge states. While for armchair edges there 
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Fig. 1: Mean (solid blue) and 
typical (dashed red) DOS for 
zigzag (left column, N z = 6) 
02 and armchair (right column, 
N a = 10) GNRs of width W = 
1.1 nm and periodic boundary 
conditions in longitudinal di- 
rection. In each panel p ty is 
given for L = 213 (1064) nm 
by red dashed (dark-red long- 
dashed) lines for 6 x 10 4 reali- 
sations of disorder. These sys- 
tem sizes correspond to 10000 
(50000) lattice sites for the 
armchair and 10392 (51960) for 
the zigzag case. For compari- 
son the mean DOS for pristine 
GNRs is given in grey in the top 
panels. 
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Fig. 2: Time evolution of 
the normalised particle density 
N\ip(ri)\ 2 for disordered GNRs. 
Device dimensions (only one 
half is shown): (1.1 x 213) nm 2 
corresponding to 6 x 1732 
atoms (zigzag) and 10 x 1000 
atoms (armchair). Times are 
measured in units of the inverse 
hopping element t = 



Fig. 3: Normalised LDOS 
(Pi/p me )U=o for particular zigzag 
GNRs with 256 x 64 sites as ob- 
tained by exact diagonalisation. In 
contrast to the other panels, where 
the LDOS is calculated at the 
unperturbed Dirac point, E — 0, 
panel (c) refers to E~ = — 0.007t. 
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is a gap at the Dirac point, E — 0, a variety of degenerate edge states exist for zigzag edges 
(see top panel of Fig. 1). 

Accounting for adatoms that lead to random on-site potentials A, a "copy" of the DOS 
centred around A will appear. The random superposition of the local A- and B-atom DOS 
results in the total (mean) DOS presented in Fig. 1. Increasing the value of A enhances the 
random on-site fluctuations, and at A = 6i (bottom panel of Fig. 1) we enter a split-band 
regime, where two distinct sub-bands emerge. 

The decay of p ty with increasing system size indicates that the states are localised in principle 
for any of the shown parameters; actually p ty vanishes for N ~ 5 x 10 4 already in all panels 
except for those in the top most row. There, however, the finite value of p ty solely indicates 
that the localisation length is larger than the system size or comparable to it. 

Figure 2 shows the time evolution of an initially localised state, as calculated by the Cheby- 
shev method. After an initial, fast spreading process (i < 10 4 io) the wave function becomes 
quasi-stationary, i.e., there are temporal amplitude fluctuations on individual sites but the over- 
all region of sites having finite amplitudes remains constant. Moderate ribbon lengths and small 
values of A result in a localisation length larger than the system size (top panel of Fig. 2) and 
therefore cause a "metallic" behaviour of the GNR. In contrast, for larger A the wave packet 
remains localised even for rather small systems. Clearly the wave function stays more localised 
for p = 0.6 < p c than for p = 0.8 > p C) since due to the lack of a spanning cluster in the 
first case, the spreading of the wave function depends heavily on quantum tunneling, whose 
efficiency decreases as A increases. 

To illustrate the localisation properties of our binary-alloy GNR model in more detail, we 
present ED results for the normalised LDOS in Fig. 3, focusing thereby on energies close to the 
(unperturbed) Dirac point, E = 0. Although for a small potential difference A and a nearly 
equal concentration of A and B atoms the wave function spans the entire (finite) lattice, the 
bimodal distribution still shows up by two interpenetrating regions with large and small wave- 
function amplitudes [see Fig. 3 (a)]. A stronger asymmetry in the concentration of type A and 
B atoms then leads to a localisation of the wave function in small regions near the GNR's 
edge [cf. Fig. 3(b)]. As may be expected the wave- function localisation effect becomes more 
pronounced at larger A, i.e. the localisation length clearly decreases upon increasing A [see 
Fig. 3 (d) and (e)]. Note, however, that for p = 0.8 a spanning cluster exists, that is we observe 
a true quantum percolation effect. 

A particularly interesting feature is the complete change in character of the eigenstates 
when crossing the energy of the unperturbed Dirac point for slightly differing on-site potentials 



and asymmetric distribution of the atom variants. While for E = + [Fig. 3(b)], the state is 
clearly localised, states on the opposite side of the impurity sub-band, E — 0~, are extended 
[Fig. 3(c)]. This transition is absent for a less pronounced asymmetry between the two atom 
types, e.g. for p = 0.6 in [Fig. 3 (a)]. 

In conclusion, we have argued theoretically and demonstrated numerically by zero-tempera- 
ture exact diagonalisation calculations that disorder of binary-alloy type in combination with 
quantum percolation effects will strongly affect electron transport in graphene nanoribbons, up 
to the point of a disorder-induced localisation of the charge carriers. We argue that in order to 
corroborate such kind of quantum localisation, transport measurements should be performed 
at much lower temperatures than used so far, but even so the effect might be covered by the 
large localisation lengths compared to the spatial dimensions of actual GNRs. 
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